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We study B s and D s spectroscopy in quenched lattice QCD using the Fermilab approach to heavy quarks. 
We obtain results at four lattice spacings, a, using 0(a)-improved Wilson quarks. We compare and contrast 
the various methods for heavy quarks on the lattice, discussing discussing which methods work best for different 
physical systems and the ease with which calculations may be performed. 



1. METHODS FOR HEAVY QUARKS 
ON THE LATTICE 

There are several ways of approximating heavy 
quarks in lattice QCD calculations, including 
Nonrelativistic QCD §]-g| (NRQCD), the static 
approximation , and the approach developed 
at Fermilab ||] which takes the uncorrected Wil- 
son action as its leading approximation but adds 
correction operators which end up resembling 
those of NRQCD rather than those of the stan- 
dard Symanzik improvement of the Wilson action 



NRQCD is based on an expansion in nonrela- 
tivistic operators (rotationally invariant but not 
Lorentz invariant) similar to that used in cal- 
culating relativistic corrections in the hydrogen 
atom. It can be thought of as arising from a 
discretization of the action arising from a Foldy- 
Wouthuysen-Tani transformation of the quark 
fields analogous to the one used in atomic physics. 



ip — > exp(9D i 'yi)i() 
leads to the action 
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The pole mass term does not affect the dynamics 
of nonrelativistic systems and is conventionally 
dropped. 



In B physics, the simplest of all the methods 
can be used, the static approximation. This ap- 
proximates the heavy quark propagator with a 
simple Wilson line moving in the time direction, 
giving an action corresponding to the first term 
in Eqn. |2|. It is clearly most useful for the heavi- 
est quarks. It is not much used recently because 
it has a much worse signal to noise ratio than 
NRQCD, which is clear in retrospect but was not 
foreseen. 

The third method can be thought of as arising 
from a partial FWT transformation: 



ip -> exp(9'Diji)ip, 
and 
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where 8' < 8, and a\ and are between and 
1. This appears to be a crazy thing to do, pro- 
ducing an action which combines the defects of 
the transformed and the untransformed actions, 
which indeed it does. On the other hand, it turns 
out that this is the action we have been stuck with 
for a long time. The Wilson action has just this 
form at leading order. The ipD 2 -^ term added to 
cure the doubling problem also contributes to the 
kinetic energy, as in Eqn. ^. Furthermore, per- 
haps surprisingly, the a\ and 02 implied by the 
Wilson action have the desired property that the 
nonrelativistic ipD 2 ip term takes over automat- 
ically from the Dirac style kinetic energy term 
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iptfhf) as k — > in the Wilson action. The Wilson 
action automatically turns into a nonrelativistic 
action in the large mass limit. 

The Wilson action does have the unwanted 
property that for large masses the pole mass does 
not equal the kinetic mass governing the energy 
momentum relation 

1 dE 

2M km ~d^' { } 

In nonrelativistic systems, the pole mass does not 
affect the dynamics and the kinetic mass governs 
the leading important term. Therefore, the pole 
mass term is normally simply omitted from the 
action in NRQCD, although there is no harm 
(and no benefit) in including it. Likewise, the 
Wilson action can be used without problems for 
nonrelativistic systems as long as M kin and not 
the pole mass is used to set the quark mass. It 
is easy to find a Wilson style action which does 
satisfy M po i e = Mki n by letting the hopping pa- 
rameter for the time direction, Kt , differ from the 
one in the spatial directions, k s . 

The Wilson action corrected up to 0(a) and 
the action of NRQCD both use one-hop time 
derivatives. When the quarks are heavy, this is 
a requirement, since two-hop time derivatives in- 
troduce new ghost states and imaginary energies 
when ma > 1. Therefore, the Wilson action ap- 
plied to heavy paths cannot follow the conven- 
tional Symanzik program of using two-hop cor- 
rections at 0(a 2 ), but must follow the NRQCD 
style of corrections which correct only the spatial 
interactions. The existence of the Hamiltonian 
and the transfer matrix ensures that this is pos- 
sible. 

The parameters of the action in this approach 
must have nontrivial mass dependence for large 
masses, just as do those of NRQCD. For large 
masses, the wave function normalization, the re- 
lation between the mass and the hopping pa- 
rameter, etc., are completely different from their 
ui — values. When m < 1 (in lattice units), this 
mass dependence may be expanded in power se- 
ries. For the Wilson and Sheikholeslami-Wohlert 
actions this yields just the usual series of op- 
erators, with the same coefficients: 

C = mipip + aim 2 ipip + . . . 



+ tplftp + bxmipflip + b 2 m 2 ip]ptp + . .. 

+ Cswipa^F^ip + C2mipcr IJ , v F^ip + ... (6) 

Thus, while for m > 1 the action becomes very 
similar to NRQCD in its behavior, for m < 1 it 
may be regarded as an all orders in m resumma- 
tion of the usual operators of the Wilson and SW 
actions. For hadrons containing charmed quarks, 
it is possible to do calculations with the Wilson 
and SW actions even with the old interpretation 
of the coefficients. However, since am 5aAQCD> 
the ability to sum up the required series in m ex- 
actly is likely to produce a much faster approach 
to the continuum limit. 

2. HEAVY LIGHT MESONS 

Here we use this approach to calculate the spec- 
trum of B s annd D s mesons. The strong interac- 
tions of a meson containing a single heavy quark 
simplify in the infinite mass limit. The heavy- 
quark spin and the light-quark total angular mo- 
mentum, j, are separately conserved. States with 
total spin J = j ± 1/2 form degenerate-mass dou- 
blets. Hence, the hyperfine splitting between the 
vector, J = 1, and the pseudoscalar, J = 0, is 
zero. Heavy quark symmetry is only approximate 
at finite quark mass so the hyperfine splittings for 
D and B mesons are proportional to I/wiq. 

The sensitivity of the hyperfine splitting to the 
heavy-quark mass makes spectroscopy of the B 
and D mesons an important check on the proce- 
dure used to determine the bare charm and bot- 
tom masses. The bare masses are important in- 
puts to our programs to determine the renormal- 
ized quark masses and weak matrix elements of 
the D and B mesons P,^0|| ■ 

Heavy-quark masses in this study are deter- 
mined in quarkonia. The bare mass is found by 
demanding the experimental spin-averaged me- 
son mass match the kinetic mass determined 
from the lattice energy-momentum dispersion re- 
lation Quark masses may also be deter- 
mined by studying heavy- light mesons. These 
masses will differ from the masses in -onia since 
lattice spacing errors differ for the two meth- 
ods |Q. Mass errors with either method how- 
ever lead to bounded uncertainties in quantities 
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Table 1 



Quenched gauge configuration ensembles. 



beta 


6.1 


5.9 


5.7 


5.5 


# configurations 


100 


350 


300 


500 


c sw = (plaq)~ 3/4 


1.46 


1.50 


1.57 


1.69 


o _1 (1P-1S) (GeV) 


2.56 (tit) 


1.81(7) 


1.16(3) 


0.75(H) 


length (fm) 


1.9 


1.7 


2.0 


2.1 



such as hadron mass splittings and matrix el- 
ements. Moreover, quark-mass uncertainties in 
these quantities are overcome by extrapolation to 
zero lattice spacing. 

3. CALCULATION DETAILS 

We study pseudoscalar and vector mesons on 
the four ensembles of quenched Wilson glue de- 
scribed in Table |l|. Lattice spacings are de- 
termined using the charmonium IP-IS splitting. 
Spacings range from a — 0.26 to a = 0.080 fm. 
The coarsest lattice only provide a check upon 
cutoff effects; results are not used for continuum 
extrapolations. 

We use Shcikholeslami Wohlcrt (SW) quarks 
with tree-level tadpole improvement. The clover 
coefficient, c sw , is determined from the average 
plaquette (see Table |l|). Calculations are done 
directly at the charm and bottom quark masses 
using. Results are fully 0(a) improved. 

Bare strange-quark masses are determined 
from the light pseudoscalar spectrum using 
leading-order chiral perturbation theory and the 
experimental it and K masses. Meson masses 
are extracted from minimum x 2 multistate fits 
to a matrix of meson correlators having Coulomb- 
gauge IS- and 2S-smeared sources and sinks. The 
full covariance matrix is retained in fits. 

4. SPECTROSCOPY 

Fig. Q shows the B s and D s hyperfine splittings 
for the lattices in Table [l]. Results from the three 
finest lattices are extrapolated linearly to a = 0. 
The D s slope has an error bound consistent with 
zero. Hence errors from the charm mass determi- 
nation are barely discernible with present statis- 
tics. The B s slope is more significant indicating 
larger errors in the bottom quark mass determi- 
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Figure 1. The D s and B s hyperfine splittings vs. 
the lattice spacing. Linear extrapolations (solid) 
with 68% confidence bounds (dashed) are shown. 
Experimental splittings (bursts) are indicated at 
a = 0. 



nation. Reducing the lattice spacing decreases 
the B s splitting. This is consistent with an under- 
estimation of the bottom quark mass in -onia |]l2| . 

The extrapolated B s and D s hyperfine split- 
tings are consistent with experimental values. 
The central values however lie below experiment. 
This may be an indication that the quenched 
B s and D s hyperfine splittings are smaller than 
experiment. Note that the quenched hyperfine 
splittings in quarkonia lie significantly below ex- 
periment. Higher statistics are necessary to de- 
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termine if there is a similar quenching effect for 
heavy-light mesons. 
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Figure 2. B s splittings divided by D s splittings. 
Hyperfine (hfs) and spin-averaged 2S-1S splitting 
ratios are shown. The experimental hfs ratio 
(burst) is shown at a = 0. Linear extrapolations 
with errors are shown. 



The ratio of the B s hyperfine splitting to the 
D s splitting, Rhfs, is proportional to m c /rn,b ac- 
cording to heavy quark symmetry. The ratio is 
expected to be less sensitive to quenching errors 
than the individual splittings. The linear ex- 
trapolation in Fig. | yields R hfs = 0.32 ± 0.08. 
This compares well with the experimental value 
0.327 ± 0.018. This is evidence that our proce- 
dure for determining the bare quark masses leads 
to the correct continuum result. 

We also study 2S states. Excited state statisti- 
cal errors are an order of magnitude larger than 
ground state errors. Excited states are also sub- 
ject to larger systematic errors. Combined er- 
rors are comparable to the spin splittings. Hence 
we report only the spin-averaged 2S splitting, 



Figure 3. The D s and B s splittings in quenched 
QCD after extrapolation to zero lattice spacing 
(octagons). Experimental IS levels (solid) are 
indicated. A potential model estimate of the 
D s spin-averaged 2S-1S splitting (dashed) is also 
shown ]l3fl . 



m(2S-lS). 

According to heavy quark symmetry splitting 
m(2S-lS) is constant up to 0(1/ mq) corrections. 
We check this in Fig. |^ by computing the ratio 
of 2S-1S splittings for B s and D s . The deviation 
from unity is consistent with 0(1 /mq) effects. 

We summarize our continuum B s and D s spec- 
troscopy results in Fig. ||. Results are plotted in 
relation to 1/mq to emphasize the behavior of 
the splittings in the heavy quark limit. We com- 
pare to experimental values for the IS states. A 
potential model inspired estimate for m(2S-lS) is 
also shown Our lattice results are consistent 
with experiment. 

Fig. ^ shows our most recent charmonium spec- 
trum, which we will use in our discussion in the 
next section. It is mainly a statistical updating 
of our previous results, except for the Xc states. 
These have been done using nonrelativistic opera- 
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Figure 4. Charmonium spectrum from the (3 = 
6.1 lattice. 



tors for the first time rather than quark antiquark 
bilinears. This allows us to obtain a result for the 
Xc 2 f° r the first time. 

5. WHICH METHOD IS BEST WHERE? 

The correction operators of NRQCD are simple 
and easy to organize. NRQCD quark propagators 
can be calculated almost instantaneously. Heavy 
quark actions which start from the Wilson action 
have the advantage that they continue to work 
as the a — ► limit is taken. However, correcting 
them to 0(a 2 ) or to (D(v A ) is messier than for 
NRQCD. 

The action of NRQCD is an expansion around 
to = co and is better and better behaved, the 
larger the quark mass. It breaks down in two dis- 
tinct ways as 1/to becomes larger. When l/(ma) 
gets larger than 1, perturbative calculability of 
coefficient functions breaks down due to factors 
of l/(ma) in the diagrams. When Aqcu/to gets 
larger than 1, the convergence of the series of non- 
relativistic operators breaks down. This series 



is very convergent for the b quark, and NRQCD 
works very well in B and T physics. The tp system 
is reasonably nonrelativistic, too, as shown by 
the success of potential models. NRQCD works 
reasonably well for this system, too, although in 
some cases not as well as originally hoped. 

For example, a recent investigation of the 
effects of v e corrections in NRQCD studied their 
effect on the hyperfine splitting of charmonium. 
In results of a few years ago, the v A , a 2 improved 
NRQCD results lay much closer to the experi- 
mental answer than the 0(a) improved relativis- 
tic Fermilab results. The interpretation seemed 
to be that the higher order corrections included 
in the NRQCD calculations made these results 
much more accurate than the less improved rel- 
ativistic results. However, last year Trottier cal- 
culated some of the v e terms which contribute 
to the hyperfine splitting and found results which 
made a big jump back down toward zero. The ac- 
tual situation seems to be that the v 2 expansion 
is very slow to converge for some quantities in the 
charmonium system. The ability of the harder to 
improve relativistic formalism to recover a Dirac- 
like form of the action in the a — > limit is an 
advantage here. 

For b quark physics, the situation is some- 
what reversed. Here, the nonrelativistic expan- 
sion is well-converged after a few orders. (Poten- 
tial models estimate that v 2 is about 1/4 in the 
J/ip and about 1/10 in the T.) On the other 
hand, these corrections have to be included to 
achieve high accuracy. One does not have the op- 
tion of replacing them and letting the relativistic 
iplpip terms in the action of Ref. || take over 
when ma < 1. This would require an absurdly 
high lattice spacing. 

A case where this plays a more significant role 
than had been expected is in the quantity M-f — 
2Mb- p5| | In the quarkonium systems, the quark 
momentum p — > mi,a s as nib — * oo. Therefore, 
(pa) 2 errors and (p/m) 2 errors also go to infinity 
in absolute terms, but approach fixed percentages 
of leading order physical quantities. In heavy- 
light systems on the other hand, the b quark's 
momentum p ~ ^-qcd as mi, — > oo. (pa) 2 errors 
and (p/m) 2 errors due to the b quark approach 
constants in absolute terms. 



6 



Of the pole and kinetic masses in the equation 
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Mkin is more affected by v 2 errors than M po i e 
since it is higher order. For example, we can ob- 
tain the spin-averaged 2Mhl — Mhh splittings 
in the D s — T and B s — ip systems from the data 
from Figs. 1-4. We obtain 0.81 GeV in the 
charm system, compared with 1.08 GeV exper- 
imentally. We obtain -1.68 GeV in the bottom 
system compared with 1.35 GeV experimentally. 
Kronfeld used potential models to show that this 
was roughly what was to be expected. jl^] The 
analogous effect in NRQCD improved only to 
0{a) is about 15%. |§ 

It should be stressed that as a percentage er- 
ror, no error in B or T physics is blowing up as 
m — > oo. For example, the uncertainty induced in 
determinations of Jb from the 20% ambiguity of 
using Mb or Mr to determine the b quark mass 
is only of order a few per cent, as can be seen by 
examining a graph of /s^/m vs. 1/rra. It is sim- 
ply one of the p 2 errors known to be present in 
the 0(a) improved theory. On the other hand, it 
is probably reduced from p T to p 2 B by obtaining 
nib from Mb rather than My in this case. It has 
been argued many times that quarkonia are nice 
systems to use for error analysis and parameter 
setting because of the tools that are available for 
understanding the errors. When those tools tell 
us that some errors in quarkonia are rather large, 
that strategy may sometimes have to be altered. 

For pure B meson physics (not considering the 
T system), the situation is probably not as bad 
for "relativistic" heavy quark actions. (Relativis- 
tic is in quotes because for b systems, relativistic 
ipflip terms in the action have almost no effect 
compared with the nonrelativistic tpT) 2 ip terms.) 
In these systems, v 2 correction terms are proba- 
bly negligible since v 2 ~ (AQci)/ m fc) 2 < 1%- It 
is therefore possible that an 0(a) improved heavy 
quark theory such as the heavy SW action may do 
an adequate job here. Perturbative corrections to 
decay constants, etc., are more annoying to cal- 
culate than for NRQCD or the m = Wilson 
theory, but by no means intractable. 
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